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Through the introduction of a class of appropriate translationally invariant trial wave functions, 
we show that the strong correlations in the lowest Landau level (LLL) reflect in finite systems the 
emergence of intrinsic point-group symmetries associated with rotations and vibrations of molecules 
formed through particle localization. This quantal molecular description is universal, being valid for 
both bosons and fermions, for both the yrast and excited states of the LLL spectra, and for both 
low and high angular momenta. "Quantum-fluid" physical pictures associated with Jastrow-type 
trial functions are shown to be reducible to the molecular description introduced in this paper. 
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Motivation. — Following the discovery [l( of the frac- 
tional quantum Hall effect (FQHE) in two-dimensional 
(2D) semiconductor heterostructures under high mag- 
netic fields (B), the description of strongly correlated 
electrons in the lowest Landau level (LLL) developed into 
a major branch of theoretical condensed matter physics 
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Most recently, the burgeoning fields of semiconductor 
quantum dots and rapidly rotating tra ppe d ultracold 

neutral gases [H E EEEl El EI EE 0, M EE E0, El 

have generated significant interest pertaining strongly 
correlated states in the lowest Landau level. Further- 
more, it is anticipated that small (and/or mesoscopic) 
assemblies of ultracold bosonic atoms will become tech- 
nically available in the near futur e I20L|21| and that they 
will provide an excellent vehicle 15|, llo , Il9l . EI El El 
for experimentally reaching exotic phases and for testing 
the rich variety of proposed LLL trial wave functions, in- 
cluding the Jastrow-Laughlin (JL) [2], composite-fermion 
(CF) [7|, Moore- Read [4j, and rotating Wigner molecule 
(REM) [1,0 trial functions. 

A universal physical description of the full LLL spectra 
(including both yrast [23[ and all excited states), how- 
ever, is still missing. To remedy this, a unified theory 
for the LLL spectra of a small number of particles valid 
for both statistics (i.e., for both bosons and fermions) is 
introduced in this paper. The LLL spectra are shown 
to be associated with fully quantal and strongly cor- 
related ro- vibrational molecular (RVM) states, i.e., with 
(analytic) trial functions describing vibrational excita- 
tions relative to the set of special yrast states that exhibit 
enhanced stability and magic angular momenta, and are 
referred to as cusp states. The cusp states, important 
as they are, represent only a small fraction of the LLL 
spectrum. The molecular trial functions associated with 
them are purely rotational (i.e., vibrationless) and were 
introduced for the case of electrons in Ref. Q under the 
name rotating electron molecules (REMs). The corre- 
sponding analytic bosonic trial functions for cusp states 
[called rotating boson molecules (RBMs)] are introduced 



in this paper; see Eq. ©. 

It is remarkable that the numerical results of the 
present theory were found in all tested cases to agree 
within machine precision with exact-diagonalization 
(EXD) results, including energies, wave functions, and 
overlaps. This numerical behavior points toward a deeper 
mathematical finding, i.e., that the RVM trial functions 
for both statistics provide a correlated basis (see below) 
that spans the translationally invariant (TI) subspace 
[H of the LLL spectrum. An uncorrelated basis, with- 
out physical meaning, built out of products of elemen- 
tary symmetric polynomials is also known to span the 
(bosonic) TI subspace [25(. We are, however, unaware of 
any other correlated functions which span this subspace. 
Indeed, although the Jastrow-Laughlin function is trans- 
lationally invariant, its quasi-hole and quasi-electron ex- 
citations are not j3(. Similarly, the compact composite- 
fermion trial functions are translationally invariant [l8l |. 
but the CF excitations which are needed to complete the 
CF basis are not [icj |. The shortcoming of the above 
well known correlated LLL theories to satisfy fundamen- 
tal symmetries of the many-body Hamiltonian represents 
an unsatisfactory state of affairs, and the present paper 
provides a remedy to this effect. At the same time a dif- 
ferent physical picture is being established: that is, the 
RVM functions and corresponding molecular point-group 
symmetries are superior in expressing and interpreting 
the emergent many-body correlations of the LLL states. 

Theory. — The RVM functions have the general form 
(within a normalization constant): 



$f iM (n 1 ,n 2 )Q%\0>, 



(1) 



where {n\,n2) indicates the molecular configuration 
(here we consider two concentric rings) of point-like 
particles with n\ in^ particles in the first (second) 
ring. The particles on each ring form regular polygons. 
The index RXM stands for either REM, i.e., a rotat- 
ing electron molecule, or RBM, i.e., a rotating boson 
molecule. $£ XM (ni, 712) alone describes pure molecu- 
lar rotations associated with magic angular momenta 
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C = Cq + niki + n 2 k 2l with k±, k 2 being nonnegative 
integers; £ = N(N — l)/2 for electrons and £q = 
for bosons. The product in Eq. (fT]) combines rotations 
with vibrational excitations, the latter being denoted by 
Q™, with A being an angular momentum; the superscript 
denotes raising to a power m. Both $^ XM and Q™ are 
homogeneous polynomials of the complex particle coordi- 
nates zi, z 2 , ■ ■ ■ i zjsr, of order C and Am, respectively. The 
total angular momentum L = C+Xm. Q™ is always sym- 
metric in these variables; <J>^ XM is antisymmetric (sym- 
metric) for fermions (bosons). |0 >= rii=i ex P[™ z i z 4 */2]; 
this product of Gaussians will be omitted henceforth. 

The analytic expressions for the <J>^ EM (for fully po- 
larized electrons) were derived in Ref. Q employing a 
two-step method: (i) First a single Slater determinant 
[that breaks the rotational (circular) symmetry] was con- 
structed using displaced Gaussians as electronic orbitals, 
i.e., 

u(z, Zj) = -= exp[-|z - Zj\ 2 /2] exp[-i(xYj + yXj)].(2) 

The phase factor is due to the gauge invariance. z = 
x — iy, and all lengths are in dimcnsionlcss units of 
IbV%, with the magnetic length Ib = y/K/(m e U) c ); 
uj c = eB/(m e c) is the cyclotron frequency. The cen- 
ters Zj = Xj + iYj, j = 1, 2, . . . , N of the Gaussians are 
the vertices of the regular polygons in the (771,772) geo- 
metric arrangement, (ii) A subsequent step of symmetry 
restoration was performed using the projection operator 

V ( C ) = h Jo" d7e l7(i ~ C) , where L = J^Zi U is the to- 
tal angular momentum operator; this yielded trial wave 
functions with good total angular momenta C @, Q . 

Analytic expressions for the $^ BM (for spinless bosons) 
can also be derived using the two-step method. Natu- 
rally, in the first step one constructs a permanent out 
of the orbitals of Eq. ([2]); one also uses the equivalence 
lu c — > 2Q between the cyclotron frequency lu c (electrons) 
and the rotational frequency Q (bosons) Here we 
present as an illustrative example the simpler case of 
N = 3 (and TV = 4 in the Appendix) bosons having a 
(0, N) one-ring molecular configuration. One has (within 
a normalization constant) 



) *2 ' 3 J > 



<&£ BM (0,3) = Yl C{hM,h) Perm^ 1 

o<ii<22<; 3 

(3) 

where the symbol "Perm" denotes a permanent with el- 
ements z\ 3 , i,j = 1,2,3; only the diagonal elements are 
shown in Eq. The coefficients were found to be: 

c(h,h,h) = (jiMj (jl p A 



cos 

l<i<j<3 



(4) 



where 1 < M < 3 denotes the number of different indices 
in the triad (h,l 2 , h) and them's are the multiplicities of 
each one of the different indices. For example, for (1,1,4), 
one has M = 2 and p\ — 2, p 2 = 1. 

The <J>^ EM expressions for electrons in a (0,N) or a 
(1,N — 1) configuration are given by Eqs. (6.2) and (6.4) 
of Ref. [9(, respectively. For electrons (1) M = N in all 
instances and (2) a product of sine terms replaces the 
sum of cosine terms appearing in Eq. ((4|) . 

We note that ^ XM (n 1 ,n 2 ) = for both bosons and 
electrons when C =/= Cq + n\k\ + n 2 k 2 . This selection rule 
follows directly from the point group symmetries of the 
(711,712) molecular configurations. 

The vibrational excitations Q\ are given by the same 
expression for both bosons and electrons, namely, by the 
symmetric polynomials: 



N 

E 



(5) 



where z c — (1/N) J2i=i z i * s the coordinate of the cen- 
ter of mass and A > 1 is a prime number. Vibrational 
excitations of a similar form, i.e., Q\ = X)t=i z i (and cer- 
tain other variations), have been used earlier to approx- 
imate part of the LLL spectra. Such earlier endeavors 
provided valuable insights, but overall they remained in- 
conclusive; e.g., for electrons in the neighborhood of the 
maximum density droplet [with fractional filling v = 1 
(u = N(N — l)/2£)], see Refs. [5| and I6J , and for bosons 
in the range < L < N, see Refs. [HI, Hi, H] . 

The advantage of Q\ is that it is translationally 
invariant (TI) j^, [r|, a property shared with both 
$^ BM and <J>g EM . In the following, we will discuss il- 
lustrative cases, which will show that the RVM func- 
tions of Eq. (fT]) provide a correlated basis (RVM ba- 
sis) that spans the TI subspace 0, H3, EH of nonspu- 
rious states in the LLL spectra. The dimension Z) TI (L) 
of the RVM-diagonalization space (using the RVM ba- 
sis) is much smaller than the dimension _D EXD (L) of 
the exact-diagonalization (EXD) [9] space spanned by 
uncorrelated determinants Det^' 1 , . . . , z£\ or perma- 
nents Perm[zi , . . . , z^} formed with Darwin- Fock or- 
bitals. The remaining D EXD (L) - D TI (L) states are 
spurious center-of-mass excitations (generated by apply- 
ing Q™) whose energies coincide with those appearing 
at all the other smaller angular momenta 0. Thus 
D TI (L) = D EXD {L) - D EXD (L - 1); see TABLE H 

Three spinless bosons — Only the (0, 3) molecular con- 
figuration and the dipolar A = 2 vibrations are at play 
(as checked numerically), i.e., the full TI spectra at any 
L are spanned by the wave functions 



(0,3)0? ^{k,m}, 



(6) 



with k, 777 = 0, 1, 2, . . ., and L — 3fc + 2m; these states are 
always orthogonal. This represents a remarkable analogy 
with the case of N = 3 electrons (see below). 
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TABLE I: LLL spectra of three spinless bosons interacting via a repulsive contact interaction g5(Zi — Zj). 2nd column: 
Dimensions of the EXD and the nonspurious TI (in parenthesis) spaces (the EXD space is spanned by uncorrelated permanents 
of Darwin- Fock orbitals). 4th to 6th columns: Matrix elements [in units of g/(TvA 2 ), A = yJh/(mQ,)\ of the contact interaction 
between the correlated RVM states {k, m} [see Eq. @]. The total angular momentum L — 3k + 2m. Last three columns: 
Energy eigenvalues from the RVM diagonalization of the associated matrix of dimension D TI (L). There is no nonspurious state 
with L = 1. The full EXD spectrum at a given L is constructed by including, in addition to the listed TI eigenvalues [D TI (L) in 
number], all the energies associated with angular momenta smaller than L. An integer in square brackets indicates the energy 
ordering in the full EXD spectrum (including both spurious and TI states). Seven decimal digits are displayed, but the energy 
eigenvalues from the RVM diagonalization agree with the corresponding EXD-TI ones within machine precision. 



L 


D EXD (D TI ) 


{k, m} 


Matrix elements 






Energy eigenvalues (RVM diag. or EXD-TI) 





1(1) 


{0,0} 


1.5000000 






1.5000000(1] 


2 


2(1) 


{0,1} 


0.7500000 






0.7500000(1] 


3 


3(1) 


{1,0} 


0.3750000 






0.3750000(1] 


4 


4(1) 


{0,2} 


0.5625000 






0.5625000(2] 


5 


5(1) 


{1,1} 


0.4687500 






0.4687500(2] 


6 


7(2) 


{2,0} 


0.0468750 


0.1482318 








{0,3} 


0.1482318 


0.4687500 




0.0000000(1] 0.5156250(4] 


7 


8(1) 


{1,2} 


0.4921875 






0.4921875(4] 


8 


10(2) 


{2,1} 


0.0937500 


0.1960922 










{0,4} 


0.1960922 


0.4101562 




0.0000000 0.5039062(6] 


12 


19(3) 


{4,0} 


7.3242187xl0 -4 


1.0863572xl0" 2 


1.5742811xl0~ 2 








{2,3} 


1.0863572x10" 2 


0.1611328 


0.2335036 








{0,6} 


1.5742811xl0~ 2 


0.2335036 


0.3383789 


0.0000000 0.0000000 0.5002441(13] 



TABLE U provides the systematics of the molecular 
description for the beginning (0 < L < 12) of the LLL 
spectrum. There are several cases when the TI subspace 
has dimension one and the exact solution coincides with 
a single {k,m} state. For L = the exact solution coin- 
cides with $R BM = 1 (Q0 = 1); this is the only case when 
an LLL state has a Gross-Pitaevskii form, i.e., it is a sin- 
gle permanent [see |0) in Eq. {!])]. For L = 2, we found 
$^ act cx Q 2 (for the index [i], see caption of TABLE HJ), 
and since [see Eq. ((SJ)] Q 2 oc (zi — z c )(z2 — z c ) + (zi — 
z c)( z 3 ~ z c) + i z 2 — z r)( z z — z c), this result agrees with 
the findings of Refs. [13j, |26j concerning ground states of 
bosons in the range < L < N. For L = 3, one finds 

$exact K $RBM_ Smce $ RBM ^ _ _ _ ^ 

[see Eq. this result agrees again with the find- 

ings of Refs. [ill, [2^|. For L = 5, the single non- 
spurious state is an excited one, <E>K?j act oc ^f BM Q 2 - 



For L = 6, the ground-state is found to be <j)?^ act oc 



-160$jP M /9- 



1/4 



fa/"* - (21-32) (zi- ^3) (22-23) , i-e-, 
the bosonic Laughlin function for v = 1/2 is equivalent 
to an RBM state that incorporates vibrational correla- 
tions. For L > N(N - 1) (i.e., v < 1/2), the EXD yrast 
energies equal zero, and with increasing L the degener- 
acy of the zero-energy states for a given L increases. It 
is important that this nontrivial behavior is reproduced 
faithfully by the present method (see TABLE QJ . 

Three electrons — Although unrecognized, the solution 
of the problem of three spin-polarized electrons in the 
LLL using molecular trial functions has been presented 
in Ref. [27] . Indeed, the wave functions in Jacobi coordi- 
nates in Eq. (18) of Ref. 27 1 are Q precisely of the form 



<5>^ M Q2 l 7 as can be checked after transforming back to 
cartesian coordinates. It is noteworthy that Laughlin did 
not present molecular trial functions for electrons with 
N > 3, or for bosons for any N. This is done in the 
present paper. 

Four electrons — For N — 4 spin-polarized electrons, 
one needs to consider two distinct molecular configura- 
tions, i.e., (0,4) and (1,3). Vibrations with A > 2 must 
also be considered. In this case the RVM states are not 
always orthogonal, and the Gram-Schmidt orthogonal- 
ization is implemented. 

Of particular interest is the L = 18 case (y — 1/3) 
which is considered [2| as the prototype of quantum- 
liquid states. However, in this case we found (see the 
Appendix) that the exact TI solutions are linear superpo- 
sitions of the following seven RVM states [involving both 
the (0,4) and (1,3) configurations]: |1 >= $f 8 EM (0,4), 



TABLE II: N = 4 LLL electrons with L = 18: Expansion 
coefficients in the RVM basis (labelled by the \i >'s) for the 
three lowest-in-energy EXD-TI states (labelled [1], [2], [4]; 
see the Appendix). The 4th column gives the RVM expansion 
coefficients of the corresponding Jastrow-Laughlin expression. 



RVM 


EXD-TI [1] 


EXD-TI [2] 


EXD-TI [4] 


JL 


|1 > 


0.9294 


-0.3430 


0.0903 


0.8403 


|2 > 


-0.1188 


-0.0693 


0.8930 


-0.1086 


3 > 


0.0067 


0.0382 


-0.2596 


0.0076 


4 > 


0.0137 


0.0191 


-0.0968 


0.0395 


5 > 


0.2540 


0.8486 


0.1519 


0.4029 


6 > 


0.0211 


0.0283 


0.3097 


0.0616 


7 > 


-0.2387 


-0.3935 


0.0877 


-0.3380 
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N=4e L=18 

[1] r ^ - P] 



Q 

X 





pure (0,4) REM 
(£=18) 



pure (1,3) REM 
(£=18) 



(0,4) REM (£=14) - 
(X =2, m=2) 



FIG. 1: CPDs for N = 4 LLL electrons with L = 18 (y = 
1/3). Top row: The three lowest-in-energy EXD-TI states 
(see the Appendix). Bottom row: The RVM trial functions 
associated with the largest expansion coefficients (underlined, 
see TABLE Hill of these three EXD-TI states in the correlated 
RVM basis. See the text for details. The solid dot denotes 
the fixed point rrj. Distances in nm. 



12 >= $* EM (0,4)Q|, |3 >- * 10 



$* EM (0,4)Q|, |4 >= 



3^(0, 4)Qi, |5 >= $? EM (1,3), |6 >= $^ M (1,3)Q|, 
and 17 >= $ REM 



15 (1,3)Q3. The expansion coefficients of 



the three lowest-in-energy EXD-TI states (labelled [1], 
[2], [4]; see the Appendix) in this RVM basis are listed 
in TABLE [TTJ One sees that for each case, one compo- 
nent (underlined) dominates this expansion; this applies 
for both the yrast state (No. [1]) and the two excita- 
tions (Nos. [2] and [4]). To further illustrate this, we 
display in Fig. [1] the conditional probability (pair corre- 
lation) distributions (CPDs) P(r, ro) (see Eq. (1.1) in 
Ref. [9|) for these three EXD-TI states (top row) and for 
the RVM functions (bottom row) corresponding to the 
dominant expansion coefficients. The similarity of the 
CPDs in each column is noticeable and demonstrates that 
the single RVM functions capture the essence of many- 
body correlations in the EXD-TI states. Full quantita- 
tive agreement (within machine precision) in total ener- 
gies can be reached by taking into consideration all seven 
RVM basis states (see the Appendix; also for other L's, 
including L — 30 with 19 RVM basis states). Naturally, 
a smaller number of RVM states yields intermediate de- 
grees of high-quality quantitative agreement. 



The celebrated JL ansatz FJ 1 



<i<j<N 



(Zi - Zj) 2p+1 has 



been given exclusively an interpretation of a quantum- 
fluid state However, since the RVM functions span 
the TI subspace, it follows that any TI trial function (in- 
cluding the JL ansatz above and the compact CF states) 
can be expanded in the RVM basis. As an example, we 
give in TABLE [ill (4th column) the RVM expansion of 
the JL state for = 4e and L = 18. One sees that, com- 
pared to the EXD yrast state (1st column), the relative 
weight of the pure (0,4) REM (denoted by |1 >) is re- 



duced, while the weights of higher-in-energy vibrational 
excitations are enhanced. In this context, the liquid char- 
acteristics are due to the stronger weight of molecular 
vibrations which diminish the rigidity of the molecule. 

Conclusions. — The many-body Hilbert space corre- 
sponding to the translationally invariant part of the LLL 
spectra of small systems (whether fermions or bosons, 
and for both low and high angular momenta) is spanned 
by the RVM trial functions introduced in Eq. ([I]). The 
yrast and excited states for both short- and long-range 
interactions can always be expressed as linear superposi- 
tions of these RVM functions. Thus the nature of strong 
correlations in the LLL reflects the emergence of intrin- 
sic point-group symmetries associated with rotations and 
vibrations of molecules formed through particle localiza- 
tion. We stress the validity of the molecular theory for 
low angular momenta, where "quantum- liquid" physical 
pictures @, 0, 03 have been thought to apply exclu- 
sively. Our analysis suggests that liquid-type pictures, 
associated with translationally invariant trial functions 
(e.g., the JL and compact CF functions), are reducible 
to a description in terms of an excited rotating/vibrating 
quantal molecule. 

Work supported by the U.S. DOE (FG05-86ER45234). 



APPENDIX 



RBM analytic expression for N = 4 bosons — Here we 
present the case of N — 4 bosons having a (0, 4) one-ring 
molecular configuration. One has (within a normaliza- 
tion constant) 



2 BM (0,4) = 



h+h+h+h=C 

£ c{h,hMM)^HA\z l 2 \z l i,z\% (7) 

0<h<l-2<l 3 <U 

where the symbol "Perm" denotes a permanent with ele- 
ments z' 3 , i,j = 1,2,3, 4; only the diagonal elements are 
shown in Eq. ([7]). The coefficients were found to be: 



M 



c(h,h,h,k) = IP** n 



Pk- 



vfc=l 



(3?i + Z 2 -^3 -3L0- 



(8) 



where V n is an operator that generates the nth permu- 
tation of the boson labels (subscripts of Z/s) 1, 2, 3, and 
4. The index M (with 1 < M < 4) denotes the number 
of different indices in the tetrad h,h, li) and the p^'s 
are the multiplicities of each one of the different indices. 
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TABLE III: LLL spectra of four spin-polarized electrons interacting via the Coulomb repulsion e 2 /{uTij) . Second column: 
Dimensions of the full EXD and the nonspurious TI (in parenthesis) spaces (the EXD space is spanned by uncorrelated deter- 
minants of Darwin-Fock orbitals). Last three columns: Energy eigenvalues [in units of e 2 /(kIb)] from the diagonalization of 
the Coulomb interaction in the TI subspace spanned by the trial functions $£g^.(0,4)Q^ and $6+3fc(M)Q™ (RVM digonal- 
ization). Third to sixth columns: the molecular configurations (n\,n2) and the quantum numbers k, A and m are indicated 
within brackets. There is no nonspurious state with L — 7. The full EXD spectrum at a given L is constructed by including, 
in addition to the listed TI energy eigenvalues [D T1 (L) in number], all the energies associated with angular momenta smaller 
than L. An integer in square brackets indicates the energy ordering in the full EXD spectrum (including both spurious and 
TI states), with [1] denoting an yrast state. Eight decimal digits are displayed, but the energy eigenvalues from the RVM 
diagonalization agree with the corresponding EXD-TI ones within machine precision. 



L 


D EXD (_D TI ) 


[(ni,Tia){fc, A, m}] 




Energy eigenvalues (RVM diag. or 


EXD-TI) 


6 


1(1) 


l(0,4){0,A,0}j 




2.22725097(1] 






8 


2(1) 


[(0,4){0,2,1}] 




2.09240211(1] 






9 


3(1) 


[(1,3){1,A,0}] 




1.93480798[1] 






10 


5(2) 


[(0,4){1,A,0}] [(0,4){0,2,2}] 




1.78508849(1] 


1.97809256(3] 




11 


6(1) 


[(1,3){1,2,1}] 




1.86157215(2] 






12 


9(3) 


[(0,4){1,2,1}] [(0,4){0,2,3}] [(1,3){2,A,0}] 


1.68518201(1] 


1.76757420(2] 


1.88068652(5] 


13 


11(2) 


[(1,3){1,2,2}] [(0,4){1,3,1}] 




1.64156849(1] 


1.79962234(5] 




14 


15(4) 


[(0,4){2,A,0}] [(0,4) {1,2,2}] 


[(0,4){0,2,4}] 


1.50065835(1] 


1.63572496(2] 


1.72910626(5] 






[(1,3){2,2,1}] 




1.79894008(8] 






15 


18(3) 


[(1,3){3,A,0}] [(1,3){2,3,1}] 


[(1,3){1,3,2}] 


1.52704695(2] 


1.62342533(3] 


1.74810279(8] 


18 


34(7) 


[(0,4){3,A,0}] [(0,4){2,2,2}] 


[(0,4){1,2,4}] 


1.30572905(1] 


1.41507954(2] 


1.43427543(4] 






[(0,4){0,2,6}] [(1,3){4,A,0}] 


[(1,3){2,2,3}] 


1.50366728(8] 


1.56527615(11] 


1.63564655(15] 






[(1,3){3,3,1}] 




1.68994048(20] 







For example, for (2,2,2,5), one has M = 2 and p\ = 3, 
P2 = 1; for (0,0,0,0), one has M = 1 and p\ = 4; for 
(1,2,3,9), one has M = 4 and p\=P2 = P3 =p± = l. 

Expressions for any number N of bosons and any 
molecular configuration (ni,n2, . . . ,n q ) can also be de- 
rived and will be presented in a future publication [2^ ]. 

LLL spectrum of four polarized electrons — For N = 4 
spin-polarized electrons, one needs to consider rovibra- 
tional states [see Eq. (1) in the main text] for two dis- 
tinct molecular configurations, i.e., $6+4^(0, 4)Q™ and 
$^(1,3)Q£\ TABLE HI] summarizes the quantal 
molecular description in the start of the LLL spectrum 
(6 < L < 15 and L = 18). 

In several cases the nonspurious (TI) states are given 
by a single trial function as defined in Eq. (1) of the main 
text. Indeed for L = 9 the yrast state is a pure REM 
state, i.e., $^ EM (1,3). For L = 11 the single nonspuri- 
ous state is the first excited state in the full spectrum, 
coinciding with the molecular vibration <&{^ EM (1, 3)Q2- 

Of particular interest is the L = 18 case; it corresponds 
to the celebrated v = 1/3 fractional filling, which is con- 
sidered [3] as the prototype of quantum liquid states. 
However, in this case we found (see TABLE [FF[| that the 
EXD-TI nonspurious solutions are linear superpositions 
of seven molecular states involving dipole (A = 2) and 
octupole (A = 3) vibrations relative to both the (0,4) 
and (1,3) configurations. Focusing on the yrast state 
with L — 18, we found that its largest component is 
the pure $^ 8 EM (0,4) REM state with a 0.9294 overlap 
with the EXD solution (see TABLE II in the main text); 
the contributions of the remaining six states are much 



smaller, but they bring the overlap to precisely unity. 
Unlike the v = 1/2 case of bosons, we stress that the 
fcrmionic Jastrow-Laughlin functions at all i/s exhibit 
less-than- unity overlaps @, Q • 



Of great interest also is the L = 30 [y = 1/5) 
case, which in the composite-fermion picture was found 
to be susceptible to a competition [29J between crys- 
talline and liquid orders. However, we found that the 
exact nonspurious states for L = 30 are actually lin- 
ear superpositions of the following 19 [— D Tl (L — 
30)] RVM functions: $^™(0, 4)g^ 2 " 2fe , with k = 

: 3 I t(l,3)g2 2_3fc/2 , with k - 2,4,6; 



0,1,2,3,4,5,6; $«f 3 ^(l, 3)Q^- 3fc/2 , v,nh /,- 



«(0,4)Q 



8 - 4fc / 3 , with k = 0,3; and $£^(1, S^"*, 



with k = 2, 3, 4, 5, 6, 7, 8. Diagonalization of the Coulomb 
interaction in the above TI subspace yielded an en- 
ergy 0.25084902 e 2 /(Kl B ) per electron for the yrast 
state; this value agrees again, within machine precision, 
with the EXD result. The most sophisticated variants 
of the composite-fermion theory [including composite- 
fermion diagonalization (CFD), composite-fermion crys- 
tal (CFC), and mixed liquid-CFC states 0, 0, HI] fall 
short in this respect. Indeed the following higher ener- 
gies were found [H [U : 0.250863(6) (CFD), 0.25094(4) 
(mixed), 0.25101(4) (CFC). The CFD basis is not 
translationally invariant [10| . Consequently, to achieve 
machine-precision accuracy, the CFD will have to be per- 
formed in the larger space of dimension D EXD (L = 30) = 
169. 
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